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ABSTRACT 

We calculate the torque on galaxies in clusters due to gravity and to dynamical friction 
forces in order to study the possible origin of small-scale alignment effects as the result 
of interactions with their environment. The equation of motion for the position angle 
of a galaxy is derived by using a pie model. We find that weak radial alignment effects 
can be produced by this mechanism involving only the most massive galaxies. We 
also introduce a dependence on the cluster eccentricity to our equations in order to 
explore the alignment of galaxies with the cluster's major axis. We find that in the 
inner regions of high eccentricity clusters, alignments of massive galaxies with the 
cluster's major axis dominate over the radial ones. This mechanism could account for 
the observed alignment effects of the most massive galaxies with the major axis of 
their host cluster. Our results suggest that dynamical friction is a viable generator of 
alignment only for the most massive cluster galaxies. For the observed alignments of 
normal galaxies a primordial origin has to be explored. 
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1 INTRODUCTION 

The observational detection of different kinds and degrees 
of non-random orientation effects of galaxies and clusters 
has a long and rich history which goes back in time to 
the end of the 19th century (see Djorgovski 1983 for a re- 
view). Their prevailing importance lies in the fact that they 
naturally arise in various cosmological models of structure 
formation, with more preeminence in models with a flatter 
power spectra on large-scales ( West et al.|1989 Dekel|19"88 1. 
In these models, galaxy formation occurs in collapsing flat- 
tened structures where colliding density shock waves will 
tend to produce spin vectors aligned with the protocluster 
plane. In fact, |Doroshkevich| ( |197d| predicted the existence 
of anisotropic orientation of galaxies as evidence for struc- 
ture formation driven by a neutrino-dominated universe. He 
showed that the average bound region around a high density 
peak in the initial Gaussian density fleld is an ellipsoid with 
principal axes parallel with those of the deformation tensor 
at the peak within the formed fllamentary structure. At the 
scale of galaxy clusters, the "peak-patch" model seems to 
validate such predictions in the context of cold dark matter 
universes ( [Bond et al. 11996] ). It is also very important to un- 
derstand how alignment effects encode the influence of the 
environment on the formation and evolution of galaxies, but 
also how they contaminate the weak lensing measurements. 



major axis of the Perseus-Pisces supercluster coincides with 
the peak of the distribution of the position angle of galax- 
ies. Djorgovski ( 1983 1 found the same effect for the Coma 



cluster and similar results have been reported for the Local 



supercluster (|MacGillivray and Dodd| 


1985| |Flin and God- 


lowski 1986 Kashikawa and Okamura 


19921. Furthermore, 



the evidence, for alignment of clusters' major axes over tens 
of megaparsecs ( [Joeveer et aLjpTSl iBinggcU 1982! [West] 
]i989j could reinforce the predictions of alignment effects for 
flatter spectra on large scales. 



On small-scales, Fong et al. ( 1990 1 found some statisti- 
cal evidence for three types of alignment: i) between neigh- 
bouring galaxies in a cluster; it) between the major-axes of 
galaxies and the cluster's major-axis; and iii) between the 
position angle of a galaxy and the radial vector to the clus- 
ter centre. In this paper we are interested in the alignments 
of galaxies within clusters such as those repo rted by |Haw- 
[ley and Peebles] ( 11975t , iThompson] ( 11976[ ) and jAdanis et al 
(19801. More recently, this effect has also been shown by 



On large scales, Gregory et al. ( 1981 1 showed that the 



Pereira and Kuhn (20051, Agustsson and Brainerd (20061, 
Torlina et al.l ( l2007| and lPereira et al.] ( ]2008[ ). They found 
statistical evidence that galaxies appear to be aligned ra- 
dially with respect to the cluster centre. Moreover, another 
point of interest concerns the evidence that in linear clus- 
ters the most massive galaxies present the same orientation 
as their parent structures ( Sastry|1968 Carter and Metcalfe 
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1980||Binggeli|1982[|Trevese et al.|1992| . One could then ask 
whether such smaU-scale ahgnments are of primordial origin 
or have arisen as the result of the environmental influence 
on galaxies residing in clusters. In the latter case, primor- 
dial tidal interactions would seem an obvious candidate to 
produce, or even to suppress any signal of anisotropy in the 



19851 



orientation (e.g. Binney and Silk|1979 Wesson|1982 Dekel 



In this paper, we use a simple analytical model to in- 
vestigate whether gravity and dynamical friction torques can 
account for stable orientations of galaxies as an evolutionary 
process. By taking into account the drag-gravity torque pro- 
duced by the central part of the cluster, we derive in section 
[2] the equation of motion for the position angle of galax- 
ies. We show that for mass galaxies of [lO^^ 10^^] Mq the 
amplitude of the oscillations decays in a Hubble time and 
discuss the viability of this mechanism to produce align- 
ment effects. A dependence of the alignment timescale on 
the cluster eccentricity is introduced in section |3] and the 
alignment of galaxies with the cluster's major axis is also 
studied. We estimate which of these kind of alignments is 
the dominant effect. A brief discussion about the implica- 
tions of our results on the theories of structure formation is 
given in section [4] 

2 RADIAL ALIGNMENTS 

In this section we quantify the effect on the dynamics of a 
galaxy, of the forces of gravity and dynamical friction, as a 
function of its distance to the cluster centre. For the cluster 



potential we adopt the triaxial model of Richstone ( 1980 1 
given by 



<f (r) = In [{x/af + {y/bf + [z/cf + l] , 



(1) 



with $c = ^TvGpc [2/3 (l/a^ + + 1/c^)] ^ , where a, 6 
and c denote the main axes length and pc is the core density 
and r = {x,y,z). The gravitational force Fg exerted by the 
cluster on a galaxy of mass M is 



F„ = -MV'I'fr) = -2Af/(r)a^ f 4* - 



where 



/(r) 



2 2 
a2 62 



+ 



-f 1 



(2) 



(3) 



and (x, y, z) are the unit directions along {x, y, z)-axis, re- 
spectively. In order to obtain a simple estimate of the rele- 
vant timescales let us assume the gravitational force due to 
a spherical cluster (6 = c = a). In this case, equation ([2| 
acquires the form 



F, = -2/(r)Mr, 
and 



fir) 



-l>c 



+ 1 



(4) 



(5) 



We shall now calculate the timescale for the position an- 
gle decay and relate it to the orbital and crossing times. We 
define the time of alignment, taiig, as the time necessary for 
the orientation of the galaxy to decay to a fraction 1/e of its 
initial value. The orientation will be specified by the angle 
formed between the major axes of the galaxy and the radius 



vector joining its centre of mass with the cluster centre, 9. A 
further assumption is that the direction of the torque acting 
on a galaxy, due to the cluster potential, is roughly constant. 
The requirements for this premise to be valid are that the 
galaxy moves in quasi-radial orbit, and that taUg ^ torb- 
Otherwise, the dynamics of the galaxy in the presence of a 
time-varying torque would be more complex, and N-body 
simulations would be required. This is beyond the scope of 
this work and will be considered elsewhere (Gonzales and 
Teodoro, in preparation). Suppose we find that taUg S> torb 
then no alignments would be observed other than those al- 
ready present in the initial conditions. However, for align- 
ments produced as an evolutionary process of the cluster, 
the number of galaxies displaying these non-random orien- 
tations must be very low under the premise that the ini- 
tial distribution of position angles is random. On the other 
hand, if it is found that taiig ^ tn, an anisotropy in the 
orientations would be easily created. Whilst for taiig ~ tn, 
there would be a diluted signal of alignments but far from 
an isotropic orientational distribution. 

In our analysis we consider interactions with the mean 
cluster potential only, and we will neglect the interaction 
with neighbouring galaxies. For this external tidal field to 
modify the orientation of galaxies a nearly homogeneous dis- 
tribution of surrounding galaxies would be required. The to- 
tal torque experienced by a galaxy is the sum of the torque 
due to the drag force and to the gravitational force. We will 
derive the expressions for the torques acting on a galaxy. 



which we model in a manner analogous to Peebles ( 1969 1 
and Wesson ( 1982 1984 1. The approximation consists in as- 



sociating a dipole to the structure of the galaxy. The dipole 
is formed by two masses m = Mgaiaxy/2, whose separa- 
tion is calculated from the position d — ^R, of the cen- 
tre of mass of half a disk (or ellipsoid) of radius (or ma- 
jor axis) R and uniform density. Thus ^ = 0.42 for disk 
galaxies and 0.37 for ellipsoidal ones. An expression for the 
torque may be calculated for galaxies when the shape is 
specified by two unequal semi-axes. The force acting on 
each mass of the equivalent dipole, according to equation 
F9''(rO ~ —2mf{vi)vi, where ri and V2 are their position 
vectors. The small variation in the value of / at the posi- 
tion of the two masses is responsible for the torque. AFg = 
-2m [/(r2)r2 — /(ri)ri]. Let 2d = r2 — ri be 



;^(2) 



the length of the dipole and 2vc 



r2 + ri the galaxy centre 



of mass galaxy. The torque is then Tg = d x AFg, where 
X denotes outer product. We obtain a useful approximation 
for the torque as a function of the angle between d and rem , 
by carrying out a Taylor expansion of AFg and retaining the 
first order term only — AFg/4m Ri /(rcm)d -I- (d ■ V/) rem- 
This leads to 



Tg = —4m (d • V/) d X Tc 



(6) 



Making use of the dipole length and of the fact that in a 
spherical potential / is a function of radius only, we obtain 



Tg = — 4mii^|/'|r- cos ^ sin ^ex. 



(7) 



where /' denotes derivative respect the radius and ex = 
ri X r2/|ri x ra]. 

In order to evaluate the dynamical friction component 
of the torque, we calculate the dynamical friction force on 
a galaxy with a typical velocity cr, where a is the cluster 
velocity dispersion. The dynamical friction expression for 
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such a galaxy can be rewritten in the form F^'^ — —kmp\, 



where k = 2n\nAG^ma~^ [erf(l) - 2e"V0r] and InA ^ 
(6 — 8.6) is the Coulomb logarithm for clusters. Thus, torque 



Td EE d X AFd 



-kmp{r)d x (v2 — vi) . 



(8) 



Note p(r) must be different for each of the masses, but since 
the galaxy size is small compared with the distance to the 
cluster centre, it is reasonable to adopt the density at the po- 
sition of the centre of mass. The velocities can be expressed 
in terms of the peculiar velocity in the cluster r and the 
angular velocity Lj = 0asvi,2 = rTdxaj which yields 



Td = —2kmpd uj. 



(9) 



Strictly, the total dynamical friction is not a sum of drag 
forces experienced for each mass: both halves of the "two- 
body" galaxy are dragged by both wakes. 

The equation of motion for the position angle is derived 
by expressing the total torque on a galaxy of characteristic 
radius R. In terms of its inertia tensor I = 2mB? and its 
angular velocity that is. 



Irl = 2mR''6. 



(10) 



If we invoke the geometric factor f = d/R, the equation 
describing the motion of the position angle of the galaxy is 
obtained from expressions 0, (|9| and ([lOf: 



6 + kp{r)fe + 2\f '{r)\r('^ sine cose ^ 0. 



(11) 



A damped harmonic oscillator analogy is evident: the major 
axis of the galaxy can either rotate or oscillate. It rotates if 
it has sufficient energy _Ecrit = [O^ — |/'(r)|r^^ cosS] /2, to 
carry it past the vertical upright position, but oscillates oth- 
erwise. If the galaxy rotates, its angular velocity is minimum 
near the position perpendicular to the radial alignment and 
maximum when is parallel to it. As the motion is damped 
by the action of the frictional term, rotation can eventually 
be replaced by oscillations around the equilibrium point cor- 
responding to the radial alignment. 

The Chandrasekahr's formula deduced under the as- 
sumption of interaction of point masses, only slightly over- 
estimates the real friction that acts on an extended body 
( Binney and Tremaine|[l987 Read et al.||2006" i. Hence, the 
simple dipole model adopted in this study should be satis- 
factory, and a description in terms of the exact geometry of 
the galaxies would not produce an important change in the 
damping timescale. The time required to damp the initial 
amplitude of the oscillations is contained in the parameters 



l3 = kp{r)C and 7 = 2|/'(r)|r^^ 



(12) 



In order to visualize whether the dynamical friction pro- 
duced in the central regions of an inhomogeneous cluster is 
strong enough to damp the oscillations in a Hubble time, we 
have compared the parameters /3 and 7^''^ as a function of 
the radius. The upper panel of figure ([TJ shows such a com- 
parison, the assumed parameters are: InA = 7, ^ = 0.4 and 
a = 0.3 Mpc. 

The gravitational term dominates in the outer regions 
of the cluster, where galaxies oscillate around the point of 
a perfect radial alignment. As the galaxy reaches the inter- 
mediate zones, the dynamical friction becomes of the same 
order of magnitude as the gravity, tending to damp the os- 
cillations. In the innermost regions of the cluster, near 600 
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Figure 1. a) Dynamical friction damping term compared with 
the gravitational term as a function of the galaxy mass, the cluster 
velocity dispersion and the distance to the cluster core; b) same 
as a) but here the gravitational term for a = 500 km/s changes 
due to its dependence on the cluster eccentricity e. In the lower 
panel the Dynamical Friction term assumes M = 10^^ ^Q- 



kpc (~ 2a) away from the centre, the dynamical friction has 
an increasing effectiveness. At this stage the oscillations are 
strongly damped, in a fraction of a Hubble timescale. 

One conclusion derived from these estimates is that any 
alignment produced by this mechanism should be very weak, 
as a result of taUg ~ tn ~ tort- Moreover, because the 
adopted values of the parameters were those associated with 
the galaxies that experienced a stronger damp, any radial 
alignment generated should preferentially involve the most 
massive galaxies in the cluster. This result is not surpris- 
ing because the dynamical friction grows in proportion to 
the mass. The observational evidence of radial alignments 
however, shows no preference for the most massive galaxies 
( |Thompson||1976[ ). A viable way of distinguishing between 
a primordial or an evolutionary origin of alignments, is by 
looking at the mass and position of the galaxies display- 
ing this effect. Alignments arising from an evolutionary pro- 
cess, should be detected in the intermediate and innermost 
regions of the clusters and involve only the most massive 
galaxies. If, at early times, the cluster had an isotropic dis- 
tribution of orientations then, after a Hubble time the dis- 
persion of the new distribution would be a factor away 
from the perfect radial alignment, showing some statistical 
evidence of alignments. Furthermore, a search for alignment 
effects by using the external galaxies could dilute a posi- 
tive signal because the isotropy in the orientations has not 
been significantly affected by the dynamical friction. On the 
other hand, if it is found that the effect extends to the outer 
cluster regions involving low mass galaxies, it will be quite 
suggestive of a possible primordial origin. 
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3 ALIGNMENT WITH THE CLUSTER 
MAJOR AXIS 

The study of this type of ahgnment is performed more con- 
veniently in cartesian coordinates. The torque on the galaxy 
is Tg = (r2 — ri) X AF5/2, where the differential gravita- 
tional force across the galaxy is obtained from equation ([2| 

AFg = -2m/(r)a — — — ^+ — ^ — .(13) 

As in the previous section, because the difference in the den- 
sity of the cluster on the position of the masses is not consid- 
erable, one can assume that the value of / is approximately 
equal to the density on the centre of mass of the galaxy i.e. 
/(r) /(ri) /(r2). Furthermore, in order to obtain some 
insight on the general feature of the motion let us restrict to 
the study the problem in the plane x — y and assume that 
b = c. Hence, the z component of the torque is 

Tg,z = -m/(r)a^ {(2:2 - xi) (y2 - yi) {b~^ - a~^) } . (14) 

Defining the eccentricity of the cluster as = 1 — {b/af , 
and 4> to be the angle between the galaxy major axis and the 
cluster's major axis, we can express the gravitational torque 
as Tg,z — —ma?f{v)(f'e^b~^ sin (/)Cos After the geometric 
factor is introduced by using the inertia tensor I = 2mR^ 
and the dynamical friction added, we obtain 

4> + kp{r)f<j) + 2/(r)C^e^(l - e^)~^ sin </> cos </> = 0. (15) 

Note that this term vanishes for e = because the alignment 
of galaxies with the "major axes" of a spherical cluster loses 
its meaning. The bottom panel of figure ([T]) shows the com- 
parison of the relative size of the parameters 

K=kp{r)f and v = 2f{r)fe'^{l - , (16) 

for different cluster eccentricities. For massive galaxies in 
clusters of low velocity dispersion, the drag and gravity con- 
tributions are approximately of the same order of magnitude 
for the inner regions. The relevance of the term containing 
the cluster eccentricity is that it projects the importance 
of the gravitational interaction. Tidal interactions between 
more eccentric shapes are stronger than those experienced 
for symmetric structures. Thus, the gravitational term domi- 
nates the oscillations of the galaxies major axes for the flat- 
test clusters. However, the observed ellipticity of clusters 
(typically e ~ 0.5) lies within the interval [0 — 0.6] where the 
dynamical friction is comparable with the gravitational ef- 
fect. One can infer then, that if such a mechanism operates in 
clusters of galaxies, in principle alignments other than those 
between the most massive galaxies with the major axes of 
the flattest clusters should be observed. Because this mech- 
anism does not damp the position angle oscillations for all 
the mass range of galaxies, one can conclude that alignment 
effects of this kind would hardly be due to an evolutionary 
process like the one considered in the present paper. In an 
asymmetric cluster there is an interplay between the radial 
and parallel alignment. A measure of which of these two ef- 
fects dominates the final orientation of galaxies located in 
the innermost regions of clusters is achieved by comparing 
the parameters 7 and v. From figure ^ we observe that 
in the outer regions of the cluster, the gravitational term 
responsible for the radial alignment is larger than that cor- 
responding to the alignment with the cluster's major axis. 
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Figure 2. Ratio of 7^/^ and w^/^: the gravitational factor respon- 
sible for the radial alignment to that responsible for the parallel 
alignment for different cluster eccentricities. The asymmetry of 
the cluster potential dominates in the innermost region of the 
cluster (~ 2a). I 

In the inner regions this situation changes and the parallel 
alignment becomes the dominant one. In this region of the 
cluster, the dynamical friction grows to a size comparable 
with the gravitational term. Therefore, the galaxy tends to 
move toward the cluster centre and to get aligned with its 
major axis. For nearly symmetric clusters, the radial align- 
ment always dominates over the parallel alignment but the 
dynamical friction is not significant, so that no alignments 
are present, except for galaxies of 10^^ — 1Q^^M0. 

4 DISCUSSION AND CONCLUSIONS 

We have studied the possibility of obtaining non-random 
orientations of galaxies in clusters produced by the com- 
bination of gravitational tidal interactions with the cluster 
potential and dynamical friction. The orientation of non- 
rotating galaxies is unstable in the presence solely of grav- 
itational tidal interactions, i.e., they oscillate around the 
position of equilibrium corresponding to the alignment be- 
tween the main axes of the galaxy with the radial vector 
to the cluster centre, or parallel to the cluster's major axis. 
With the inclusion of the dynamical friction we found that 
these oscillations can be significantly damped only for galax- 
ies with masses in excess of lO^^Mg. The damping occurs m 
a Hubble timescale, which is of the same order of magnitude 
as the orbit decay timescale. Thus, the final orientation of 
the galaxies that have time to reach the cluster centre will 
be determined by the eccentricity of the cluster. Whilst at 
the outer regions of nearly symmetric clusters the gravi- 
tational term responsible for the radial alignment is larger 
than that corresponding to the alignment with the clusters's 
major axis, in the inner regions the parallel alignment be- 
comes the dominant one where the anisotropy of the tidal 
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field is maximum. For clusters of low eccentricity, e ^ 0.6 
(ellipticity ^ 0.3), parallel alignment is only important very 
close to the core. For flatter clusters, where the tidal field is 
markedly anisotropic, the gravitational term (v) is more im- 
portant than it is in the symmetric case (7). The anisotropy 
of the tidal field is of the same order of magnitude of the ra- 
dial action of gravity even at larger distances of the cluster 
centre (see Figure[T|. In these regions, the dynamical fric- 
tion damps the oscillations around the cluster's major axis 
in a fraction of a Hubble time. These results makes plausi- 
ble that massive galaxies that have had time to reach the 
cluster centre are aligned with the cluster's major axis. One 
can then conjecture that the final result of these encoun- 
ters would produce a giant, central galaxy aligned with the 
cluster's major axis. Even when tidal disruption or merger 
with other galaxies may occur, the direction defined by the 
clusters' major axes would be preferential to the motion of 
the remnants of these processes. The tidal field defines a 
preferred direction, along the major axis for eccentric clus- 
ters, and along the radius vector to the cluster centre if it is 
nearly spherical. This mechanism could produce anisotropic 
mergers of objects with quasi-parallel vectors in the cen- 
tral region. Moreover, because of the low number of ob- 
jects of very massive galaxies in a cluster, mergers would 
be highly anisotropic and the final product would have a 
non-neglegible peculiar velocity. Such central galaxies have 



been reported by Hill et al. ( 1988 1 



There is some numerical evidence that the orientation 
of the first-ranked galaxies formed by the merging of sub- 
structure in clusters of galaxies, is determined by the ini- 
tial ellipticity of the cluster ( Rhee and Roos|1990 l and that 
possibly orientational memory is preserved by the galaxies. 
However, such a model shares the drawback of the merging 
theories which cannot account for the high peculiar velocities 
of cD galaxies. Further, the merging theory cannot account 
for alignment of galaxies with the clusters major axes in the 
external regions. On the other hand, with the mechanism of 
drag and gravity torques studied here, alignment of galax- 
ies could not be produced at the outer regions of clusters 
where the damping timescale is several times the Hubble 
time. Thus, the preference for the most massive galaxies 
at the central region of clusters suggest that gravitational 
tides combined with dynamical friction are not the only vi- 
able generator of the observed alignments, which include low 
mass galaxy alignments extending to the outer region of the 
clusters. Radial and parallel alignment effects must therefore 
be compatible with the idea of a primordial origin. However, 
one of the reasons why this effect is not strongly observed in 
all the clusters is that there are many other environmental 
effects which could suppress the alignments at a post cluster 
formation era. The final orientation of galaxies would be dif- 
ficult to determine for galaxies with non-radial orbits. These 
are some of the factors not contemplated in this work. Con- 
cerning a primordial origin, [West et al.| |T989 ) have shown 
that alignment effects of clusters of galaxies naturally arise 
in Universes dominated by Mixed Dark Matter. In fact the 
small dispersion in cD's luminosities, for a number of possi- 
ble cluster physical conditions suggest that cD progenitors 



should be special ab initi o (|Westnl994[ ). Barnes and Efs- 
|tathiou| ( |1987[ ) and |Bond| ( |1987[ ) have also suggested that 
non-random orientation of density perturbations could be 
primordial, and be present at the initial conditions. Lam- 



bas, Groth and Peebles have found evidence of alignment of 
elliptical galaxies with the major axis of their host clusters, 
but not for spiral galaxies. If the orientation of galaxies in- 
deed refiect the primordial orientation of maxima of the den- 
sity field, then the angular momentum of galaxies need not 
have any connection with the angular momentum of ellipti- 
cal galaxies, which have a lower specific angular momentum. 
These ideas are being investigated (Gonzalez & Teodoro, in 
preparation). 
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